Tomography of correlation functions for ultracold atoms via time-of-fiight images 
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We propose to utilize density distributions from a series of time-of-fiight images of an expanding 
cloud to reconstruct single-particle correlation functions of trapped ultra-cold atoms. In particular, 
we show how this technique can be used to detect off-diagonal correlations of atoms in a quasi-one- 
dimensional trap, where both real- and momentum- space correlations are extracted at a quantitative 
level. The feasibility of this method is analyzed with specific examples, taking into account finite 
temporal and spatial resolutions in experiments. 
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I. INTRODUCTION 

Ultracold atomic gas provides a controllable platform 
to study strongly correlated many-body physics, which 
has attracted strong interest recently To reveal 

many-body properties of the underlying systems, nor- 
mally one needs to detect some kinds of correlation func- 
tions. The detection method for atomic systems is typ- 
ically very different from detection of condensed mat- 
ter materials. For solid state systems, linear responses 
provide a useful method to measure the system corre- 
lation functions. For ultracold atomic gases, measuring 
linear response is possible but not always convenient. A 
powerful and widely used detection method for atomic 
gases is based on the time-of-flight (TOF) imaging tech- 
nique [H, which is unique to these atomic systems and 
has no counterpart in condensed matter materials. Dur- 
ing the TOF imaging, one measures the light absorption 
of an expanding atomic cloud released from the exter- 
nal trap. The light absorption gives information of the 
density distribution of the expanded cloud, which, under 
some approximation, is proportional to the initial mo- 
mentum distribution of the atomic gas prior to expan- 
sion. Therefore, the TOF imaging provides us a useful 
technique to extract the diagonal single-particle correla- 
tion function in the momentum space. With only the 
diagonal correlation in the momentum space, in general 
it is inadequate to reconstruct the real space correlation 
function. To fulfill this gap, several methods have been 
proposed to introduce additional control techniques, and 
measure the real space correlation functions based on the 
atomic interference 0] or the Fourier sampling Q . Some 
of the real-space correlations can also be inferred from 
measurement of the two-particle correlations which can 
be extracted from noise spectroscopy 

In this manuscript, we discuss a method to extract in- 
formation about full single-particle correlations in both 
the momentum and the real spaces, by measuring the 
density profiles in the TOF images at different expansion 
times. Our detection proposal is in the spirit of tomog- 
raphy in the spatial-temporal space, where the full cor- 



relations are reconstructed through certain algorithms. 
The method is strict for (quasi) one-dimensional (ID) 
systems, and can be applied to higher-dimensional cases 
where the symmetry is applied to reduce the effective 
dimensionality. Compared with other detection schemes, 
this method is more direct in the sense that it relies solely 
on the TOF images and does not require introduction 
of other challenging control techniques. Besides, it only 
utilizes the mean density distribution of the resulting im- 
ages, for which the detection is in general significantly 
easier than the measurement of quantum noise of the cor- 
responding images. As an example of applications of this 
detection method, we consider a quasi-lD Bose gas and 
demonstrate that all off-diagonal single-particle correla- 
tions can be reconstructed from a series of TOF images at 
various expansion times. The feasibility of this method is 
analyzed under realistic experimental conditions of spa- 
tial and temporal resolutions. 



The remainder of this manuscript is organized as fol- 
lows. In Section mi we discuss the main reconstruction 
formalism by analyzing the ballistic expansion process of 
an atomic gas and showing how the single-particle corre- 
lations are obtained from density profiles in the TOF im- 
ages. We consider both the quasi-lD case and the high- 
dimensional cases where one can either separate variables 
or have spherical symmetry to reduce dimensionality. As 
a by-product, we also derive in this section the formula 
for the far-field limit, where the initial momentum distri- 
bution is directly connected with the final density profile, 
and obtain the quantitative conditions under which this 
far-field-limit formula is valid. After introducing the gen- 
eral formalism, to illustrate its applications we give an 
example in Section IIIII by considering a quasi- ID Bosc 
gas, and demonstrate how correlations are reconstructed 
from the TOF images. Finally, we discuss the effects 
from finite spatial and temporal resolutions and conclude 
that this method is applicable within present experimen- 
tal conditions. 
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II. FORMALISM FOR RECONSTRUCTION OF 
CORRELATION FUNCTIONS 

In most experiments, the atomic gas is usually pre- 
pared inside some external optical or magnetic trap. 
When the system achieves its thermal equilibrium, one 
detects the properties of its underlying many-body state 
through the TOF imaging. To perform this TOF imag- 
ing measurement, one turns off the external trap such 
that the atomic gas starts to expand in space. For sim- 
plicity, we assume a ballistic expansion process where 
the atomic interaction is negligible during the expansion. 
This is typically a good approximation for optical lattice 
experiments when the atomic density is not high || . For 
strongly interacting atoms near a Feshbach resonance, to 
get the ballistic expansion, one needs to first sweep the 
magnetic field to the deep BEC or BCS side to turn off 
the interaction, as done in many experiments. 

For a ballistic expansion, the atomic cloud freely ex- 
pands and the expansion dynamics is already known. Be- 
fore turn-off of the trap, the system is in thermal equi- 
librium, and its single-particle correlation function is de- 
noted by G (r, r') = (^(^^(r')), where <fo(r) = 4>{r,t < 
0) is the atomic field operator, which is either bosonic or 
fermionic corresponding to bosons or fermions, respec- 
tively. Assume at time t = the trap is turned off, and 
we measure the density profile (n(r, t)) = (cj)^ (r,t)(f>(r,t)} 
of the atomic cloud at various expansion time t. From 
the measured density (n(r, t)) @, we would like to recon- 
struct the full single-particle correlation Qo(r,r'). Note 
that this reconstruction is in general impossible for a 
three-dimensional system, where Go(r, r') depends on six 
variables and (n(r, t)) has only four variables. The corre- 
lation Qo(r,r') therefore contains more information than 
the density profile (n(r, t)). However, we will show in the 
following that for any quasi-lD systems, the reconstruc- 
tion of the correlation functions can be done exactly. Fur- 
thermore, the correlation can also be obtained for many 
practical high-dimensional systems, where the system ei- 
ther possesses spherical symmetry or can be separated in 
variables for the single-particle correlations. 

The expansion dynamics of the atomic cloud is well de- 
scribed by the Schrodinger equation for the atomic field 
operator (taking H = 1): 

idt<Kr,t) =-^(r,i). (1) 
Ira 

The solution of this equation in the momentum space is 
simply given by 

0(k, t) = Wo(i)0 o (k) = o (k)e- 4fc2 * /2m , (2) 

where </>(k, t) is the Fourier transform of the field operator 
4>{v, t) at time t, b( is the free propagator, m is the atomic 
mass, and initial conditions are 0(k, t = 0) = 0o(k). The 
expectation value of density distribution thus takes the 



form 

(»(',*)) = // dkidk 2 e- i(kl - ka) - r 

xe^- k ^ t / 2m (4(k 1 )0 o (k 2 )). (3) 

It is useful to define new variables k + = ki + k2 and 
k = ki — k2, with which the density expectation value 
becomes 

WM)> - JJ ^-^"^ 

where t = t/2m is defined to simplify the notation. 

An important feature of the density profile in Eq. Q is 
that for long evolution time t (called the far-field- limit), 
the exponential term gives a rapid-oscillating phase fac- 
tor except for the region around the point (k_ = 0, k+ = 
r/t). As a consequence, the integration over k + and k 
is dominated by the contribution from such region, lead- 
ing to an approximating form for the density expectation 
value (see the derivation in Appendix A) 

Notice that the final density profile is directly propor- 
tional to the initial momentum distribution via a scal- 
ing relation k = r/(2t). This formula has been widely 
used for interpretation to the measurement result from 
the TOF images. For this far-field-limit to be valid, the 
expansion time needs to be sufficiently long. Quantita- 
tively, it has to satisfy the following condition 

(Ak+ • Ak_)t> 1. (6) 

Here, Ak+ and Ak_ are characteristic scales for the ex- 
tension (variation) of single-particle correlation function 
C/o(ki,k2) along the k+ and k_ directions, respectively. 
The derivation of this condition for the far-field limit is 
shown in Appendix [X] 

The extraction of momentum distribution only utilizes 
the measured density profile in the far-field limit. The 
density distribution measured at other expansion times 
contain more information about the single-particle cor- 
relations. Next, we discuss how to inverse the relation 
in Eq. ([3]), and to reconstruct the off-diagonal correla- 
tions from TOF images. We first discuss in detail the ID 
system where the reconstruction algorithm is exact, and 
then generalize the method to higher dimensional cases 
by reducing to a set of ID problems when symmetry or 
separability arguments are applicable. 

Let us consider a quasi-ID atomic cloud along the x- 
direction where the transverse degrees of freedom are 
frozen through a deep potential such as an optical lat- 
tice. The density profile {n(x 1 t)) thus depends on cor- 
relation functions through a double-integration with two 
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phase factors along the spatial and temporal directions, 
respectively. One can think about performing a double 
Fourier transform to inverse the relation. Specifically, we 
obtain 



/OO pOO 
dx / dt(n(x,t))e~ l( - px+UJt) 
-oo JO 

m I ,+ f—2muj/p~p\ f—2muj/p + p 
a s ^ 



2b| 
% 2m 



, I / dfc + i V o \ " I (7) 

4ir p \J + -2mw/p + k + ' 




for the cases with p ^ 0, which is connected with the off- 
diagonal momentum correlations and hence of particular 
interest. 

Notice that the integration over time is for t > only, 
since all the density profiles are obtained after turn-off of 
the trap. This partial Fourier transform hence generates 
the principal value (W) integral on the right hand side of 
Eq. j7]). This expression can be significantly simplified if 
the single-particle correlations are real. This condition is 
equivalent to assume a time-reversal symmetry through- 
out the system. In fact, for systems in a stationary frame, 
it is very unlikely for the time-reversal symmetry to be 
broken as long as there is no such symmetry breaking 
term in the Hamiltonian Q. Under this condition, the 
principal value integral is purely imaginary, and the off- 
diagonal momentum correlation (with fc_ ^ 0) can be 
solved as 



41/cJRe 



dt 



dxe -ik-{k+l-x) ( n (a. )t )) 



(8) 



Therefore, with the measured density profile (n(x, t)) at 
different time t, we can reconstruct all the momentum 
correlations through the formulae © and ||5J). With the 
knowledge of all momentum space correlations, the real 
space correlation function can also be obtained via a dou- 
ble Fourier transform. After substituting the expression 
([5]), the integration over two momentum indices can be 
carried out analytically, leading to 



(fh{xi)<j) {x2] 



2tt 



-Re 



dt / dx 

11 J -oo 



f 2 



-ix-(2x— x+)/4i 



n{x,t)\ 



Here, the variables x+ = X\ + X2 and x- = x\ — X2 
are defined to simplify notations. The convergence of 
the integration over time for t — ► is usually guaranteed 
by cancelation from the fast oscillation in the exponent 
e -tx_(2a-x + )/4t_ xhese two simple relations between the 
measured density profiles and the initial correlation func- 
tions in both the momentum space [Eq. ([8])] and the real 
space [Eq. ©J are the central results for this section. 



After introducing the reconstruction scheme for the ID 
case, next we consider the problem in higher dimensions. 
Since the density profile (n(r,t)) is in fact a function 
of d(spatial)+l(temporal) variables, in the most general 
case it is no longer possible to fully reproduce the single- 
particle correlation, which is a function of d x d variables. 
However, if the system has some properties which allows 
us to reduce dimensionality, a similar procedure can still 
be performed on the reduced problems. For instance, if 
there exists spherical symmetry in 3D systems, the field 
operator can be written as <j>(r,t) — <f>(r,t), where r is 
the radial distance. In this case, the free expansion pro- 
cess is governed by the Schrodinger equation in spherical 
coordinate, 



id t <f>(r,t) 



1 

2m 



d$ + -d r )4>{r,t). (10) 
r 



By defining a new operator ip(r,t) 
tion above becomes 



id t ip(r, t) 



d 2 
2m 



(r, t), the equa- 



(11) 



Notice that this equation has exactly the same structure 
as the Schrodinger equation fl} for the ID case, which 
allows us to obtain the single particle correlation function 

((/)J(r)0o(' ,/ )) = ( , 0o( r ')V ; o( r ' / ))/( rr ' / ) f rom the density dis- 
tribution n(r, t) via a similar procedure as Eq. ([5]). 

The reconstruction scheme can also be generalized to 
higher dimensional cases when the correlation functions 
are separable, i.e., 

{^{v,tW,t))= [] Ul(ri,t)<t>M,t)) . (12) 



In this case, the higher dimensional problem can be effec- 
tively reduced to a set of ID problems and then treated 
separately using the same method discussed above. 



III. CORRELATIONS IN ONE-DIMENSIONAL 
QUASI-CONDENSATES: AN EXAMPLE 

In the previous section, we have shown how to recon- 
struct the single-particle correlations from the measured 
density profile for ID systems and higher-dimensional 
cases where symmetry or separability can be used to re- 
duce the effective dimensionality. To illustrate applica- 
tions of this general formalism, in this section we consider 
a specific example and show hot to reconstruct the cor- 
relation functions for bosonic atoms in a quasi- ID trap. 
The reconstruction formula in Eqs. (jHJ) and @ are given 
as double integrals of the measured density profiles. In 
realistic experiments, however, one can take only a finite 
number of images, and each image contains a finite num- 
ber of resolvable data points. In the following discussion, 
we will take into account the finite spatial and temporal 
resolutions, and show that a good approximation to the 
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correlation function can be inferred already from about 
a dozen of images. 

In this example, we consider a quasi- ID Bose gas 
trapped through either a highly elongated cigar-shaped 
potential [8j or a deep transverse two-dimensional optical 
lattice [9(. In this quasi- ID geometry with strong radial 
confinement lo p , atoms at low temperatures are essen- 
tially frozen at the ground state of the radial harmonic 
trap. Thus the radial wave function is governed by 



where a a 



4> P {y,z) = J— g 



j, p = 1/ \jmujp characterizes the extension of the 
radial wave function and m is the atomic mass. When 
a p is much greater than the effective length scale of in- 
teratomic potential R e , the gas can be described with an 
effective ID interaction rate [Kj| 



-(y 2 + z 2 )/(2al) 



(13) 



SlD 



2a, 



(14) 



where a s denotes the s-wave atomic scattering length in 
free space. 

The static single-particle correlation function in this 
quasi-lD system has been calculated by evaluating phase 
fluctuation effects around the saddle point condensate 
densities. This process gives the axial real-space correla- 
tion function [ljj] 



g (x,x') « vW^oMe- 172 ^^, (15) 

where no(x) is the zero-temperature Thomas- Fermi den- 
sity distribution in the axial direction 
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FIG. 1: (Color online) The density profiles along the axial 
direction at different expansion times for a quasi- ID Bose gas 
released from the trap. Throughout the paper, we consider a 
gas of 87 Rb atoms with total number N = 10 4 interacts via 
s-wave scattering length a s — 5.45 nm. The gas is confined in 
a cigar-shaped potential with the trapping frequencies u) x = 
2tt x 30 Hz and uj p = 2tt x 20000 Hz. The Thomas-Fermi 
radius along the axial direction is J?tf ~ 160/im. Such a 
system can be realized in a 2D optical lattice with the lattice 
constant d = 426 nm and the lattice depth Vo — 10-Er, where 
E-r, is the recoil energy. In this plot, we consider a spatial 
resolution of Aa; = 20/im, which gives about 16 data points 
(squares) from the absorption image at the very beginning 
of t = 0. The temperature of the system is assumed to be 
T = 2Quj x . 
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and the function F s (x,x') takes the form 



F ' {x > x) = mz 



in 



{1-x/R tf )(1 + x'/Rtf) 



{l + x/R TF )(l-x'/R TF ) 



(16) 



(17) 



Here, T is the temperature, N is the total atom number, 
Ho denote s the chemic al potential at the trap center, and 
i?TF = v2aW ( muj x) i s the Thomas-Fermi cloud radius 
along the axial direction (with a small axial trapping fre- 
quency lo x ). We emphasize that we only use the correla- 
tion function in Eq. (fT5|) as an example to illustrate our 
detection method. The derivation of this correlation is 
irrelevant for our following purpose. 

We assume a ballistic expansion of the atomic gas along 
the axial direction for TOF imaging (the transverse trap 
is still on), which allows us to numerically simulate the 
expansion process with the initial correlation given by 
Eq. (|15[) . and obtain the expected density distribution 
at each measurement time. Using these resulting den- 
sity profiles, then we try to reconstruct the correspond- 
ing correlation functions using Eqs. © and ([5]). From 



the derivation in the last section, if we have infinite spa- 
tial resolution and take an infinite number of images, we 
should be able to obtain exactly the same correlation 
function as shown in Eq. (fT5|l . So the purpose here is ac- 
tually to analyze the effects from finite spatial and tempo- 
ral resolutions as is the case for realistic experiments. In 
the following discussion, we assume to take only 10 ~ 30 
images, and for each image we only know (n(x,t)) for a 
discrete set of points from the finite spatial resolution. 
In Fig. [TJ we show a typical set of results for (n(x, <)) at 
various times. Notice that we assume here a finite spatial 
resolution of Ax = 20/im with R^-p/ Ax w 8, which gives 
about 16 readable data points from the absorption image 
at the very beginning of t = 0. 

First, let us try to extract the correlation function 
in momentum space using the formula in Eq. ([5]). In 
Fig. [2l we plot the reconstructed off-diagonal correla- 
tion (0q (/c)0o ( — k)) using the simulated density profiles. 
There are three important features one can read from this 
plot. First, the reconstructed correlations are fairly close 
to the expected exact values obtained from the Fourier 
transform of the real-space correlation in Eq. (|15p . espe- 
cially when the correlation is sizable such that the error 
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FIG. 2: (Color online) The momentum-space correlation 
(0q(/c)0o( — fe)) reconstructed from the density profiles via 
Eq. JS), where a total of Mt = 30 images are taken equidis- 
tantly between t = and the maximal expansion time 
imax = 6.8 s (triangles) and t max = 13.7 s (dots). The results 
are compared with the expected exact values from a direct 
Fourier transform of Eq. (|15[1 (solid line). The oscillatory be- 
havior for large k is an artificial effect due to the sudden drop 
of the Thomas-Fermi density distribution at the cloud edge. 
Notice that the reconstruction fails for k — > (hollow points) , 
while the questionable region shrinks (arrows) with increasing 
tmax and Xmax- In this plot, we assume a spatial resolution 
of Ax = 40/xm with Rtf/Ax = 4, which gives about 8 data 
points from the absorption image at the very beginning of 
t = 0. Other parameters used here are the same as the ones 
in Fig. [1] 
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is relatively small. Second, the reconstruction fails for 
k very close to zero. This is because the TOF images 
taken here are not for an infinite time duration. For 
the reconstruction, we take the data points within the 
ranges given by x max and i max for coordinates and time, 
respectively. Since the reconstruction relies on a Fourier 
transform, these finite ranges set a limit for the resolu- 
tion in the momentum space. In fact, by increasing the 
evolving time t max , which simultaneously requires an in- 
crease of a; max since the cloud expands more, we can push 
the reconstruction further towards k = 0. Third, except 
for the region k ~ 0, the results are insensitive to the 
time split between two subsequent images. Notice that 
since we use M t = 30 images for both reconstructions in 
Fig. [21 the trial with longer evolving time has larger time 
split. The result, however, is fairly close to the other trial 
for k not very close to 0. 

Next, we analyze the correlation function in the real 
space. In Fig. [3Ja), we show the reconstructed corre- 
lations of ((f>l(x)(j)o(—x)) from the density profiles via 
Eq. ©. The most striking feature of this plot is that 
the spatial correlations can be obtained very precisely 
for x not too close to 0, using as few as M t = 10 images. 



FIG. 3: (Color online) (a) The real-space correlation 
(0o(a;)0o(— x)) reconstructed from the density profiles via 
Eq. (|9]|, where the maximal expansion time tmax = 6.8 s 
are sliced equidistantly to obtain Mt = 10 (triangles) and 
Alt = 20 (dots) images. Results are compared with the 
expected values of Eq. (|15[) (solid line), (b) Log plot of 
{<jy a {x)<j>o{— x)) around the center of the cloud at different tem- 
peratures. Exponential decay can be clearly observed in this 
plot, with the correlation length given by an exponential fit 
(solid lines). Other parameters used here are the same as the 
ones in Fig. [2] 



Even in the questionable region for x less than the spatial 
resolution Ace, some information can still be extracted 
by applying linear interpolation between image pixels, 
and the precision in that region can be significantly en- 
hanced when more images are used for reconstruction. 
In Fig. [3^b), the real-space correlations for two different 
temperatures are plotted in log scale, showing explicitly 
the exponential decay around the center of the trap with 
corresponding correlation lengths. 
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follow the same trend with variation of the resolutions. 
In Fig. [U we show the normalized RMS error in Eq. (fT5|) 
as functions of spatial [Fig. Ufa)] and temporal [Fig.fJJb)] 
resolutions, respectively. Notice that the error decreases 
rapidly by increasing resolutions in both variables, as one 
would expect. Besides, the value of the error saturates 
to a fairly small number (~ 5%) around Rtf/Ax ~ 4 
and M t ~ 20, respectively, indicating that the technique 
is quite feasible under the present technology. 



IV. CONCLUSION 

In summary, we have proposed a method to recon- 
struct the full static single-particle correlation functions 
in both the momentum and the real spaces for cold 
atomic gas by measuring the density profiles at differ- 
ent expansion times with the time-of-flight imaging. The 
method applies to quasi-lD systems and can be general- 
ized to higher dimensions when symmetry or separability 
arguments can be used to reduce the effective dimension- 
ality. As an example, we consider a quasi-lD Bose gas 
and demonstrate how real- and momentum-space corre- 
lations are reconstructed at a quantitative level. The 
feasibility of this method is analyzed by evaluating the 
reconstruction error with various spatial and temporal 
resolutions, and the result suggests that the correlations 
can be inferred with pretty good precision already with 
a dozen of images at practical spatial resolution. 
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FIG. 4: (Color online) Normalized RMS error of the recon- 
structed correlations for Go(x, —x), (a) as a function of the 
spatial resolution Aa; with a fixed number of images Mt = 20, 
and (b) as a function of the number of images Mt with a fixed 
spatial resolution Rtf/Ax = 4. The system is at tempera- 
ture T — 20uj x , with total expansion time t max = 6.8 s. Other 
parameters are the same as the ones in Fig. [T] 

Finally, we investigate how the increase of spatial and 
temporal resolutions can enhance the reconstruction pre- 
cision. To quantitatively evaluate the precision, we con- 
sider the real-space correlation as illustrated in Fig.[3ja), 
and define the root-mean-square (RMS) of the error 

Af i ,/ EiI'i[go(^,-^)-go TOF (^.-^)] 2 

A/rms -n (0)V N x 

(18) 

where Gj OF (xi, — £») is the reconstructed correlation 
from the TOF images, and the function is normalized 
to the center density no(0). Notice that this normalized 
RMS error will depend on the number of reconstruction 
sampling points N x , which is 40 as in Fig. [3Ja). But its 
value will not change much with different N x and will 
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APPENDIX A: FAR-FIELD LIMIT 

For simplicity of the notation, here we derive the far- 
field limit formula only for the ID case. The extension of 
the formula to higher dimensions is straightforward. In 
the ID case, the expectation value of the density distri- 
bution in Eq. Q takes the following form 

I = -^75 / dx J dye^f(x, y), (Al) 

where f{x,y) represents the single-particle correlation 
function. The correlation function f(x,y) spreads over 
certain ranges, with the characteristic length scales along 
the x and y directions denoted by AX and AY, respec- 
tively. The integration over dy in Eq. (|A1|1 can be per- 
formed after a Fourier transform, leading to 

1 f 

I = — / dx,^(x. s = xt) 
4tt J 

= 4^/ ds#(s/t,s), (A2) 
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where &{x, s) is the Fourier transform of f(x,y) as a 
function of y. Notice that from the Fourier transfor- 
mation theorem, the function & must spread along the 
x and s directions with the characteristic length scales 
given by AX and 1/AY, respectively. 

In the expression above, the integration over s is es- 
sentially along the line (s/i, s). For large enough time 
t, this line is almost the s-axis. Thus, the result can be 
approximated by 

J "4^/^ ( °' s)= 4^ /(0 ' 0) ' (A3) 

where the second step is given by Fourier transforming 
the function back to the {x, y) plane. The last equa- 
tion gives the far- field limit result. 



In order to make this approximation valid, the function 
&(s/t,s) and J^(0, s) should be close to each other for 
typical values of s. As the function ^(s/t,s) has char- 
acteristic length scales AX and 1/AY respectively along 
the x and s directions, the condition s/t <C AX must be 
fulfilled for s within the typical range given by 1/AY. 
Therefore, we conclude that the far-field limit given by 
Eq. (|A3|) is valid when the following condition is satisfied 

AXAYt > 1. (A4) 

The above argument can be easily generalized to higher 
dimensional cases, leading to the far-field limit condition 
as shown in Eq. (J6]). 
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